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I. INTRODUCTION 



The Kondo problem for dilute magnetic impurities in metals has attracted continuous 
interest for years. Many fascinating aspects of the problem have been revealed by various 
methods such as the renormalization group approach the local Fermi liquid theory 
and the Bethe-ansatz solution @J^. Recently a series of works by Affleck and Ludwig has 
shed a new light on our understanding of the Kondo effect P, pi]|JTT| . Their new machinery 
is a technique of two-dimensional boundary conformal field theories (CFT) |[T^. Various 
critical properties of the Kondo problem including the multichannel overscreened as well as 
the underscreened case have been studied extensively. 

The first crucial step made by Affleck in this approach was to recognize that for the ordi- 
nary Kondo problem the impurity effect is incorporated as modifying boundary conditions 
on the otherwise free conduction electrons in the low-energy effective theory From the 
CFT viewpoint, therefore, the impurity effect should be explicitly observed in the finite-size 
spectrum of the Kondo Hamiltonian. Assuming that the low-energy effective Hamiltonian 
is expressed in the Sugawara form of the current algebra Affleck and Ludwig verified the 
finite-size spectrum making use of the fusion rule hypothesis [0]. On the other hand, the 
Bethe-ansatz solutions to various models for the Kondo problem have been obtained 0,0. 
Thus the finite-size spectrum should be directly obtained starting with the microscopic 
models. It is then worth comparing the results with those predicted by boundary CFT. 

In this paper we address ourselves to this task and derive the finite-size spectrum in 
the Kondo problem from the Bethe-ansatz solution of variants of impurity models. We will 
show that the essential properties of the spectrum predicted by boundary CFT indeed agree 
with the exact results deduced by the Bethe-ansatz method. In Sec. II we begin with the 
Anderson model which is more tractable than the conventional s-d exchange model (Kondo 
model) to calculate the finite-size spectrum. We see that to the first order in 1/L with L 
being the system size the excitation spectrum coincides with that of free electrons modified 
by the phase shift arising from the impurity scattering. A similar observation was made 
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earlier for the s-d exchange model in [^]. Applying the finite-size scaling we then obtain the 
canonical (integer) exponents characteristic of the local Fermi liquid. We also argue that 
the anomalous exponents related to the X-ray absorption singularity are read off from the 
finite-size spectrum. On the other hand, the correlation effects due to the impurity go into 
the piece in the spectrum. Here the spin and charge susceptibilities for the impurity 
govern the excitation spectrum. In Sec. Ill our analysis is extended to the SU(z/) Anderson 
model (degenerate Anderson model). We shall point out that the excitation spectrum is 
written in terms of the Cartan matrix of the OSp(z/, 1) Lie superalgebra, which turns out to 
be responsible for the canonical exponents inherent in the local Fermi liquid. In Sec. IV we 
turn to the s-d exchange model in which the charge degrees of freedom of the impurity are 
completely frozen. We will see that a careful treatment of the charge degrees of freedom is 
necessary in this case. Some technical details are summarized in Appendix A. 

II. ANDERSON MODEL 

A. Bethe-ansatz solution 

In this section our purpose is to calculate the finite-size spectrum of the Anderson model. 
We first recapitulate the Bethe-ansatz solution of the model while emphasizing the charac- 
teristic properties common to the Bethe-ansatz solutions to exactly solved models for the 
Kondo problem. The Anderson model for the magnetic impurity is defined by the Hamilto- 
nian 

H = ^ekc\^Cka + VY^{c{^^d„ + 4cfc,a) + ^d'^dld^ + Ud\di^d\di, (1) 

with standard notations The model describes free conduction electrons coupled with 
correlated ci-electrons at the impurity via the resonant hybridization V . After reducing 
the model to one dimension with the use of partial waves, we linearize the spectrum of 
conduction electrons as = vk with v = 1 near the Fermi point, which is equivalent to the 
assumption of the constant density of states. The hopping operator in the kinetic term is 
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then replaced by —id/dx in the coordinate representation. These simphfications enable us to 
apply the Bethe-ansatz method to diagonalize the Hamiltonian 0J^. Upon diagonalization 
we are led to introduce two types of rapidity variables, kj and A^, for the charge and spin 
degrees of freedom, respectively. For the Hamiltonian (|1|) on a circle of circumference L the 
Bethe-ansatz equations were obtained by Wiegmann and by Kawakami and Okiji [[^ 

M 

k,L + 5{k,) = - Y: ei{B{k,) - Ap), (2) 

/3=1 

M N-2M 

27rJ, + ^^2(Aa-A/3)+ E Oi{A^- B{kj)) = -2Rek{Aa)L-2Re5ik{Aa)), (3) 

13=1 j=l 

where On = 2tan~-'^(2x/?7,) and 6{k) = — 2 tan~^(l^^/(A; — e^)) is the bare phase shift due to 
the potential scattering by the impurity at U = 0, which gives rise to the resonance width 
r = V'^/2. The total number of (up-spin) electrons is denoted as (M). Here k{Aa) is the 
complex solution with ImA; > of the charge rapidity which satisfies iB{k) — iA — UT = 
where 

B{k) = k{k-U -2ed). (4) 

Comparing Eqs.(0) and with the ordinary Bethe-ansatz equations we observe additional 
phase shift terms 6{k) proportional to 1/L whose existence is indeed quite relevant to de- 
scribe the physical properties of the impurity part. 

B. Ground-state energy 

The exact ground-state energy of the Anderson model is given by 

N-2M M 

^0= E k, + 2j2x{A^), (5) 

j=l a=l 



where 

a;(A) = e, + U/2 - f A + (e, + U/2y + [(A + (e, + U/2ff + UTyAf'] (6) 
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Note that this expression is divergent because of the constant density of states with charge 
rapidities kj distributed over [— cxd,/5] in the thermodynamic hmit. In order to evaluate the 
finite-size corrections to the ground-state energy, therefore, we have to regularize the sum 
(^) upon taking the large-L limit while keeping the 1/L scaling term which is expected to 
be universal. Most naively one may cut off the fc-integration; k G [—D,(3] with D being the 
bandwidth. This prescription, however, is not suitable for our purpose since introducing the 
explicit momentum cutoff in this way would break conformal invariance. 

Instead we are led to introduce a smooth cutoff function f{k) which is unity near the 
Fermi point and behaves as |A;|~''^"''^'* (e > 0) when \k\ oo. Eq.(|^) is now replaced by 

N~2M M 

Eo= J2 <^(%)A;,+2^y.(x(A,))x(A«). (7) 

j=l a=l 

We apply the Euler-Maclaurin summation formula to Eqs.(^, and to evaluate the 
finite-size corrections. The details of computations are relegated to Appendix A. The final 
result is given by 

En = LEo — , (8) 

" " 12L 12L ' ^ ^ 

where Vc and Vg are the velocities of massless charge and spin excitations. It is shown that 
these velocities take the same value v which is equal to unity in our convention (see Appendix 



A). To compare the result with the finite-size scaling law in CFT |jT3| one has to replace 
L with 21 since L has been defined as the periodic length of the system. Then we find the 
scaling behavior of boundary CFT [0 for the charge and spin sectors, respectively, but 
with the same Virasoro central charge c = 1. 



C. Excitation energy 

Our next task is to calculate the finite-size corrections to the excitation energies. For this 
it is sufficient to take the thermodynamic limit of the Bethe-ansatz solution. Introducing 
density functions p{k) and cr(A) for the charge and spin rapidities, respectively, we write the 
Bethe-ansatz equations in the thermodynamic limit. Standard calculations yield [^1^,3 
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Pi.k) = — + yA(A;) + B\k) r dXa,{B{k) - A)a(A), (9) 
a(A) + r dA'a2(A - AV(A') + f dka^{K - B{k))p{k) = A{A) + yZ{A), (10) 

Ja J—oo J-i 



where 



an{x) = — (11) 



A{A) = / dkai{B{k) - A), (12) 

J—oo 

/oo 
dkai{B{k) - A)A{k), (13) 
-oo 

and A{k) = {T/7r)[{k - eaf + T^]"!. It is clear from Eqs.(D, (|10D that the density functions 
can be decomposed as a = + {1/L)aimp and p = ph + {1/L)pijnp corresponding to the 
host and impurity contributions. For later convenience let us rewrite Eq. ( p!OD in such a way 
that the /c-integral is taken over the region oo] instead of [— oo,/5]. This can be easily 
performed by the Fourier transform. Then Eq.(|TO|) turns out to be 

poo poo 

a(A)+/ dAa2(A-A)a(A) = / dka^iA - B{k))p{k). (14) 

Ja J 13 

The number of conduction electrons is expressed in terms of the density functions of 
the host part 

M, / pocj poo \ 

-^=Ar-^ dkph{k)-2j^ dAa^(A)j, (15) 

where J\f = J^^dkph{k) is an irrelevant infinite constant. Similarly for the number of 
down-spin electrons Mh we get 

= / dAa,(A). (16) 

-L Ja 

The total energy has contributions from both of the host and impurity part 



E = L I dk 



2^ + 1^"=) 



<;.(*), (17) 



where we have introduced the dressed energies es{k) and ec{A) satisfying 

POO 

Esik) = -^{k)k + dAai(A-5(fc))£e(A), (18) 

J a 

EciA) = 2td + U+ r dkB'{k)ai{B{k) - A)es{k) - H dA'a2(A - A>e(A'). (19) 

J/3 Ja 

Let us now evaluate the finite-size corrections in the excitation energy spectrum. We 
first expand the total energy (|l^) to the second order in the variations Aa and A/5 of the 
integration limits a and (3. The first order terms give a chemical potential which is usually 
set to zero. However, since the shift of the chemical potential due to the impurity is of order 
1/L we must incorporate it as the finite-size corrections to the energy spectrum. In order 
to estimate this, it is crucial to notice the relation 
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^ -(2e, + f/) = / dAa,(A), (20) 

which holds for the host density function. This relation implies that the impurity level 
will change if we vary the number of conduction electrons by tuning a since our convention 
is to measure the impurity level from the Fermi level. Using Eqs. (|15|) and (pOD one obtains 
the energy shift due to the impurity 

where (= J^aid^da)) is the average number of localized electrons given by 

nd = l- r dAa,„p(A), (22) 



at the impurity site. Using the Friedel sum rule = 25f/t^ with dp being the phase shift 
at the Fermi level we now find 

Ai^d) = -^AAT,. (23) 



The second order terms are obtained from Eqs. (|17|) ~ (|19|). The result reads 



— cr(a)- 



da 



(A,8) 
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(24) 

k=p ^ 



It is convenient for our purpose to express A/5 and Act in terms of the change of the number 
of down-spin electrons and the magnetization Sh (= Nh/2 — M^) of host electrons. Taking 
the derivatives of Eqs.(|15|), (p!6| ) with respect to a and (3 we have 



(25) 



L da 



dk— = ah{a)^sc{(y), 



13 



da 



(26) 



L 8(3 



dA 



dp 



-PhmcsiP), 



(27) 



L da 



roo 

-ahia) + / dA 

J a 



da 



-o"/i(tt)^cc(a), 



(28) 



where the 2x2 dressed charge matrix (^j^) is given by the solutions to the integral equations 



^ssik) = 1 + / dAai(A - B{km,,{A), 



roo roo 

e.e(A) = / dkB {k)a,{B{k) - A)^Uk) - / dA'a2(A - A )U{A ), 

J 8 Ja 



(29) 



(30) 



ics{k) 



dAa,{B{k)-K)U{K), 



(31) 



roo roo 

eec(A) = 1 + / dkB {k)a,{B{k) - A)Uk) - / dA'a2(A - A )^,,{A ). (32) 

Consequently Aa and A/? are written in terms of ASh and AM^ which represent the devi- 
ations from the ground-state values. From Eqs.(|2^ ) ^ (|28|) we thus obtain 

^^(2) _ 27rt;c 1 pe^Agfe + U^M^ y ^ 2nvs 1 f 2(,,ASh + U^M^ y ^^^^ 



L 2\ dei{iij} J ^ L 2\ det{^ij} 
where the "velocities" of spin and charge excitations are defined as 



2TTVr 



1 



+ 



1 o-imp(a) 



ah{a) L ah{ay 



de^iA) 



da 



(34) 
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2tiv. 



1 1 PimpiP) 

For vanishing magnetic fields, we can solve Eqs. 



desjk) 
d(3 



k=l3 



(35) 



) ~ (^2]) using the standard Wiener- 
Hopf method. Note first that the dressed charges in these equations become independent of 
a in the limit of zero magnetic field for which j3 — > — oo. Hence we can simply set a = — oo, 
which corresponds to the so-called symmetric model with 2ed + U = 0. The dressed charge 
matrix is thus cast in a simple formula 



^ss 
^sc 



J 



V 



T.J 



(36) 



with P - 
with the kernel given by 



oo. Here ^(/3) is obtained from the solution to the equation 
1 



1 + — / dkB {k)K(B{k) - B{k ))^{k 
zn J/3 



(37) 



K(x) 



dp- 



-p/2 



COS ■ 



px 



2UT Jo '^coshp/2 2UT 

One can readily solve this equation by the Wiener-Hopf technique in the limit of /3 
obtaining 



(38) 



-oo, 



^sc 



73/ 



(39) 



It is remarkable that the dressed charges obtained here are independent of parameters such as 
the Coulomb repulsion and the electron filling. As we shall see momentarily this fact ensures 
the canonical exponents for conduction electrons inherent in the (local) Fermi liquid. 

The particle-hole type excitations are readily taken into account. Let us define non- 
negative integers n+ and to specify their contributions in the charge and spin sectors. 
Substituting now the explicit values of the dressed charge into Eq.(^) we write down the 
finite-size spectrum in the following form 



27rVr 
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2ttv, 



(40) 



Furthermore we notice that the "velocities" f c, f s can be written in a transparent form with 
the use of susceptibihties. Recall that the host and impurity parts of the charge and spin 
susceptibilities are given by 



X, 



imp 



imp 
As 



dud 

dH 
dH 



' imp ' 



a 



-4a, 
-PkiP) 

Pimp 



da 



A=a 



V da 



A=a. 
1 



dp 
dssik) 



k=f3 



dp 



k=l3 



(^kpimp{k) . 



where the magnetization of the impurity part reads 

1 

^ J — oo 

We find 



V. = V 



imp- 
(.c 



1 + 



\ -yArnp 



(41) 
(42) 
(43) 
(44) 

(45) 

(46) 
(47) 



where the velocity of free electrons is f = 1 in the present formulation. 

Combining Eqs.(p3D and (^Oj) we finally arrive at the simple expression for the finite-size 
spectrum 



E = Eq + a^(i) + ae(2) = + + ^^2, 
where the terms of order 1/L and are given respectively by 



2'nV . ^ x2 A 



and 



L2 



2-Kv r(AA^ft}2 



2'KV X 



imp 



X's 



-[{AS.r + nt]. 



(48) 



(49) 



(50) 



We remark that the bulk electrons do not contribute to the piece because of the 

linearized dispersion. Hence the terms of order 1/L and order give the finite-size 

spectrum of bulk electrons and impurity electrons, respectively. 
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D. Critical exponents and universal relations 



We now wish to discuss our results (^)~(^) in view of the finite-size scahng in CFT. It 
has already been found that the low-energy critical behavior in the Kondo effect is described 



by boundary CFT in which we have only left (or right) moving sector of CFT ^10,111. The 
1/L formula (^) in fact exhibits the scaling behavior predicted by boundary CFT. Notice 
again that when comparing with CFT formula one has to replace L with 21 as pointed out 
before. We observe clearly from Eqs.(|^) and (^) that the charge sector is described by 
c = 1 CFT with U(l) symmetry and the spin sector by c = 1 SU(2) Kac-Moody CFT at 
low energies. 

In Eq.(|49|) it is seen that the finite-size spectrum of the bulk electron part depends on 
the non- universal phase shift 5f- We should remember that this constant phase shift 5f 
is equivalent to the chemical potential due to the impurity, see Eq. (|2^) . This imphes that 
the effect of the phase shift amounts to merely imposing twisted boundary conditions on 
conduction electrons P, p!0| , pJ]| . When reading off from Eq. (^9D the dimension of the scaling 
operator associated with conduction electrons, we should thus ignore the dp dependence. 
Then the scaling dimension x of the conduction electron field is determined by choosing the 
quantum numbers 

/^Nh = 1, A^;, = 1/2, (51) 

which yields x = 1/2 upon ignoring Sp- This is consistent with the fact that the system is 
described by the strong-coupling fixed point of the local Fermi liquid 

According to boundary CFT |12|, however, there exists a one-to-one correspondence 



between energy eigenstates and boundary scaling operators. Such operators have scaling di- 
mensions explicitly dependent of the phase shift. Our intention now is to identify a boundary 
operator which is specified in particular by a "single-electron" quantum number (0). This 
type of excitation could be relevant to describe the long-time asymptotic behavior of the 
electron correlator in which the hybridization V in Eq.(|ll) is switched on suddenly at time 
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t = and V{t) = V for t > 0. The correlation exponent will be obtained by setting ANh = 1 
and ASh = 1/2 in Eg. (^9]) , but without ignoring 6f- We get 

ris = l- 25f/t^. (52) 

Actually such phenomena have been well known to occur as the orthogonal catastrophe 
which is induced by the time-dependent local perturbation, e.g. a power-law singularity in 
the X-ray absorption profile [Q. It is quite remarkable that the critical exponent for the 
excitonic correlation part is precisely given by Eq. (^) . Therefore we observe interesting 
evidence that boundary CFT will also play a role to clarify the critical behavior related to 
the orthogonal catastrophe. 

We point out another important property characteristic of the local Fermi liquid, i.e. the 
universal relations among the static quantities. To see this we first recall that in Tomonaga- 
Luttinger liquids (c = 1 CFT) there exist universal relations among the velocities and the 
susceptibilities such as |16|,OJl^,|19| 



VcXc = , VsXs = — , (53) 

71 /vr 

for zero magnetic fields, where Kp is a Gaussian coupling constant which features the U{1) 
conformal critical line in the charge sector. An important point is that the above relations 
hold in the present case not only for the conduction electron part but also for the impurity 
part with the fixed value Kp = 1 independent of the interaction. We can see this from the 
fact that the term proportional to 1/L^ of the energy spectrum (^) with n'^ = = can 
be cast into the expression analogous to the 1/L term 



^(2) _ o,..»mp(^^rf) 



2 



AEl^lp = 2vrt;r^^-^ + 2nvrP{As,y. (54) 



^imp — - c ^ 

Here we have introduced the velocities of the impurity part through 



imp _ ^ imp _ ^ 

„ imp' „ imp-' \^^J 

O'Jc O'Js 



with 7*"^^, 7*™P being the specific heat coefficients of the charge and spin degrees of freedom 
of the impurity electrons. This expression is essentially the same as in the bulk. From (^2|), 
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(^) and (|54D it is shown that the universal relations for static quantities (^3]) with Kp = 1 
hold for the impurity part. This is nothing but the properties of the local Fermi liquid 0. 



III. SU(zy) ANDERSON MODEL 

An integrable generalization of the infinite-f/ Anderson model with SU(z/) spin degrees 
of freedom is given by the Hamiltonian 0] 




+vjzj dx(5(x)(|m)(0|c^(x) + 4(x)|0)(m|), (56) 

m=l 

where impurity electrons have z/-fold degeneracy. |0) and \m) (m = 1,2, ■■■,//) represent 
the unoccupied and single occupied states of the impurity site, respectively, and the double 
occupancy is forbidden by infinite repulsive Coulomb interactions among impurity electrons. 
The Hamiltonian ([56|) is presented in the coordinate space and the spectrum for conduction 
electrons has already been linearized. This model has often been used to investigate the 
Kondo effect for rare-earth impurities in a metal. 

The diagonalization of the SU(z/) Anderson model can be performed by introducing 
u — 1 spin rapidities A^'^ (/ = 0, 1, ■ ■ ■ , z/ — 2) besides the charge rapidity A^'^^^^ According 
to Schlottmann |p the Bethe-ansatz equations take the form 

(/ + 1)AWl = 2vrJ« - [vr + 2tan-i((A« - e^)/(/ + 1)1)] (57) 

U-l Plq Kg 

+2 E E E tan-i((A« - Af)/{1 + q - 2p)r), 

g=0p=0/3=l 

a = 1,2,- --.Ki, Ki = Ni-Ni+i, = 0, / = 0, 1, ■ ■ ■ , z/ - 1, 

where integers /^'^ are quantum numbers associated with the spin and charge degrees of 
freedom, and Ni is the number of electrons with each spin. 

We henceforth concentrate on the excitation spectrum. As mentioned in the previous 
SU(2) case, we can work with the rapidity density functions cr^')(A) defined in the thermody- 
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namic limit when calculating the finite-size corrections to excitation energies. These density 
functions satisfy i/-coupled linear equations [§] 

,(iHA)-^ + ^ I ^ + (581 

^ ^" 27r +7rL(A-e;)2 + [(/ + l)rp ^^^^ 

I^— 1 Pig „^ 

EE / 'dAV^)(A')ir«(A-A'), 



q=0 p=0 ' 



/ = 0,l,---,i/-l, 



where piq = min(Z, q) — 6iq and F = The integral kernel is defined by 

^1"?^^^ ;rx2 + [(/ + g-2j9)r]2- 

We note here again that the density functions are decomposed into the sum of the host part 
and the impurity part: cr'-'-' = cr^''' + (l/L)(Tj-,2p- The total number of host electrons is then 



given by 



L 



E(^ + l) rdAa«(A). (60) 

1=0 -^-"^ 



On the basis of these equations we compute the finite-size spectrum following the method 
explained in the preceding sections. To proceed with clarity we will treat the charge and 
spin sectors separately. 



A. Charge excitation 

Let us first evaluate the finite-size corrections in the charge sector. For this purpose we 
can put Bi —>■ —OO for / = 0,1, — 2 while B,y^i = Q, which corresponds to the case of 
zero magnetic fields. The value of Q is determined by the number of host electrons through 
Eq. (|5DD . The total energy including the impurity contribution is simply expressed as |Q 

L J-oo [2n ttL (A - e/)2 + (z/r)2j ^ ^' ^ ^ 

in terms of the dressed energies e^^\A) (/ = 0, 1, ■ ■ ■ , z/ — 1) which are defined by the following 
equations 
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.«(A) = uK5,,., f dA'irir'HA - A')e^'-'\A'). (62) 

The finite-size corrections are computed by expanding Eg . (1611) in AQ to the second order. 
Similarly to the SU(2) model, the first order term in AQ arises from the impurity part. 
Expanding the energy to the linear term in AQ we obtain 

^E^'^ = -i^^Q = -^-Ag, (63) 

OQ IT 

where we have used the Friedel sum rule 

dE 

2i^6f = 2tt— = 27mf, (64) 

with the average number n/ of locahzed electrons given by 

% = ^ r Aat-'\A). (65) 

We next rewrite AQ in terms of the change of the host electron number AN^. From 
Eg . (|60D one finds 

l~QQ=''^h W) + '^y_^dA — =a}^ '{QKu-i,u-i{Q), (66) 

where the dressed charge ^u~-i,u~i is determined by 

e.-i,.-i(A) = ^ - r dA' E <-l^(A - A')e.-i,._i(A'). (67) 



The Wiener- Hopf solution to this eguation yields ^y^i^y^i{Q) = ^ 0. Hence we get 

\Q)L 

Furthermore the host electron density function a'^^ at A = Q is similarly evaluated as 
^'h ^^^^Qi) — \fvl2'K. We now obtain the corrections linear in 1/L 

Ai?(i) = -^AAT,. (69) 

Turning to the second order terms we expand the energy as 
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dQ 



mi 

A=Q 2 



(70) 



Using (|6^) we obtain 



L 2v 



(71) 



where the velocity is defined by 



2'KVr = — 



1 



^(-1) 



+ 



1 a. 



{u-l) 
imp 



dQ 



(72) 



Combining Eqs. (|7TD and (^), and repeating similar manipulations as in the SU(2) case we 
obtain the finite-size spectrum for the charge sector 



E = Eq + —El + -r^E^ , 



(73) 



where the terms of order 1/L and are given respectively by 



L ^ L 



2u 



nv ( 5f 
— v\ — 



(74) 



and 



2tiv xT^ \{,^Nhf 



2v 



(75) 



L2 

Here the velocity of free electrons is f = 1 and a non-negative integer represents the 
particle- hole excitation in the charge sector. Putting = we can express the 1 / L^-order 
term in an analogous fashion to the 1 / L-order term 



L2 



2z/ 



(76) 



where the velocity of the impurity part is f, 



imp 



h I imp 



B. Spin excitation 

In order to deal with the spin sector we take the integer- valence limit, nj ^ 1 {Q ^ oo), 
to facilitate our calculations. In this limit the Bethe-ansatz equations for the spin part 
reduce to [SI 
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a«(A) 



u-2 



27T (A+l/J)2 + l/4 



J2 dAA,,(A-AVW(A') 



q=0 -^^1 



/ = 0, 



(77) 



where J = — T/ej and ^^^(A) is the Fourier transform of 



\x\/2 



(78) 



The values of i?;'s are determined by the total number of electrons Ni p. We have 



L L 



/ dAa^^-^\A)-2 dA(T(''(A)+/ dAa(^+^)(A) 



(79) 



for / = 1, 2, ■ ■ ■ , z/ — 3 and 



^-^ = l-2rdA.I»)(A)+rdA.<'>(A), 

L L JBo JBi 



-2 K-1 _ _2 



L L 

We define the dressed energy function 



dAa^''-^\A) + / dAa('^-3)(A). 

B^-2 JB^-3 



(80) 



(81) 



.(0 



SO as to express the total energy in the form 



iyASi,,^2 - E r dA'Ai,{A - A')e'-'^\A'] 

n ■J Bn 



(82) 



E 



oo 1 

dA- 



1 



£(°)(A). 



(83) 



L Jbo 27r(A + l/J)2 + 1/4 
We observe here that Eqs.([77|) ~(^) take the form common to integrable SU(z/) symmetric 



models such as the SU(z/) Heisenberg model P0| , pi| and the degenerate Hubbard model |^ 
Based on the analysis of these SU(z/) models let us introduce the following quantities 



M, 



(0 



'dAa«(A) 



/ = 0,1,- ■■,1^-2, 



(84) 



L JBi 

as quantum numbers characterizing the spin degrees of freedom. Making use of the SU(i/) 



results 21,22] we find the finite-size spectrum of the spin sector 



E — En 



2n 



v-2 



1=0 



(85) 
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where AM^ = (AM/°^ 



h 5 



AM, 



(u-2). 



V = diag(fo, f 1, ■ ■ ■ , f,y_2) and nf' are non-negative 



integers representing particle-hole type excitations in the spin sector. The dressed charge 
matrix is defined as {Z,y)im = ^im{Bm) where 

u-2 



6™(A) = 5im-J2 dA'A,^(A - A')6,(A'). (86) 
We note that the matrix Zi, enjoys a nice property of {Z~^)^ Z~^ being equal to the Cartan 



matrix for the SU(z/) Lie algebra [^. In Eg . (|85D the velocities of spin excitations have been 
defined by 



Vl 



1 l^i^Bi)] de('\A) 



+ 



OB, 



A=Bi 



/ = 0,l,---,j/-2, 



(87) 



which are independent of the spin index / for zero magnetic fields. Consequently the finite- 
size spectrum for the spin sector reads 



27rvl 



2tivT^ 1 



E-E^ = ^^AMlCAMh + — |— -(LAm^)C(LAm), 



L 2 ""'^ ' ' L2 2 

where we have set n/" = for simplicity, C,, is the (z/ — 1) x (z/ — 1) SU(z/) Cartan matrix 

/2 -1 Q\ 

-1 ••. ■•. 

-1 

-12/ 
and Am denote the variations of the (z/ — l)-component vector m defined by 



(88) 



v 



(89) 



m«= £^dAa£^(A), 



(90) 



for the impurity part. This completes the calculation of the spin excitation spectrum. 



C. Finite-size spectrum 

We are now ready to present the final result for the finite-size spectrum. The 1/L piece 
in the spectrum is obtained by combining the charge and spin excitations given by Eqs.([71D 
and (|88|) . Setting ?t,+ = n'^ = for simplicity we write the result in the u x u matrix form 
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-El 



AM^C/AM-^i.(^) , (91) 



L L 2 

where AM^ = (AM^, AM^''"^^) with AM^ = AA^;, - uSp/n. Here the u x u matrix Cf 
is given by 

-1 ■-. ■-. 

C/ = . (92) 

■•• 2 -1 

U -1 ij 

It is amusing that this matrix coincides with the Cartan matrix for the OSp(z/, 1) Lie su- 
peralgebra. This matrix characterizes the energy spectrum of free electron systems as well 
as the electron models based on the OSp(i/, 1) superalgebra. That is, the OSp(i/, 1) Cartan 
matrix in the finite-size spectrum gives rise to canonical exponents for the correlation func- 
tions. In fact the critical exponents for electron correlators (c](t)cj(0)) ~ t~^" are obtained 
by specifying the quantum numbers AM^ = 1 (or 0) for Z > n (or I < n). We then find 
the critical exponents r]n = 1 irrespective of n (= 0, 1, ■ ■ ■ , z/ — 1). We see that the finite-size 
spectrum for the host part (1/1/ order) has the form in accordance with boundary CFT with 
U(l) symmetry in the charge sector and level-1 SU(i/) Kac-Moody symmetry in the spin 
sector. 

On the other hand, the 1/L^ part does not take such a simple form since the spin and 
charge velocities are different from each other due to the electron correlation effects. It is 
worth, however, considering the hypothetical situation where these two velocities are equal, 
say v*™^ = v*"'^ = 1. Then the excitation energy is again cast into the form of the OSp(i/, 1) 
Cartan matrix 

A^;^^) = §i(LAm^)C/(LAm). (93) 

Since the critical exponents of correlation functions do not depend on the velocities this 
result indicates that the spectrum of scahng dimensions are solely characterized by the 
Cartan matrix. Furthermore note that the spectrum for the impurity part {^/L'^ order) is 
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essentially given by the same formula as that for free electrons except for the modification 
of velocities. This result nicely fits with the local Fermi-liquid picture of Nozieres 0. 



IV. S-D EXCHANGE MODEL 

This section is devoted to the analysis of the finite-size spectrum of the traditional s- 
d exchange model. The reader may be somewhat curious about why the analysis of the 
simplest model comes at the end. As we shall proceed it is seen that the simplest nature of 
the model gives rise to peculiarity when trying to extract the finite-size spectrum based on 
the Bethe-ansatz solution. 

The Hamiltonian of the s-d exchange model is given by 

H = J:J dxcUx){-^-^)c^{x) + JJ:S■ci{0)cT^^,c^,{0), (94) 

where conduction electrons are scattered by the impurity spin located at the origin and the 
exchange coupling is assumed to be antiferromagnetic. Applying the Bethe-ansatz technique 
Andrei and Wiegmann [^] obtained the Bethe-ansatz equations 



M 

k,L = 27ilj + J2 ^(2A/3 - 2), J = 1, 2, ■ ■ ■ , iV^, (95) 

M 

Nhe{2K - 2) + ^(2A„) = -27rJ, + e{K - A^), « = 1, 2, ■ ■ ■ , M, (96) 

(3=1 

where e{x) = -2tan-i(x/c) with c = 2J/(1 - J^), and Nh, M are the number of host 
electrons and the number of down spins, respectively. In contrast to the Anderson model 
the charge degrees of freedom for the impurity are completely frozen, so we need a careful 
treatment of the charge excitation in the Bethe-ansatz equations. A peculiar aspect in the 
Bethe equations is that the charge and spin degrees of freedom are completely decoupled. 
Hence a standard technique of the dressed charge cannot be applied in a straightforward 
manner to the charge sector of the model. For example, a naive calculation of the dressed 
charge yields ^ = 1 for the charge sector, which is the value for spinless fermions. This seems 
apparently in contradiction to ^ = for free electrons. 
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Thus we deal with the charge degrees of freedom directly without relying on the dressed- 



charge approach. The total energy E is derived from the Bethe-ansatz solution p3|,|2l 
When there are no magnetic fields we have 

E = -rT.^j+Nh dA[^^(2A-2)-7r]a(A), (97) 

L -^-^ J-oo 

where the density function a (A) obeys 

"'^> - 2^{t (2A - 2)^ + + I (2AFT^ - L (A - AT + c^ "'^ 'I' 

with B = —oo for zero magnetic fields. 

Let us first examine the charge excitation. In the limit B = — oo one can solve Eq.(p8D 
using the Fourier transform. Having obtained the explicit form of (t{A) it is not difficult to 
evaluate the 1/L contribution in the charge sector. One finds 

1 . ^ 27CV 1 , . ^r, 71V , ^ 

_^Et = —-(AN, - If - -. (99) 

Notice that the shift 6Nh SNh — 1 originates from the tc/2 phase shift. 

Consider next the spin excitation spectrum with Nh being kept fixed. The spin part of 
the total energy in the presence of a magnetic field is [ P^ , ^ 



POO 

E' = nJ [e{2A - 2) - 7r]a(A), (100) 
Jb 

where B is determined through the number of down spins of host electrons which is 
given by 



Jh 
L Jb 



dAa/,(A). (101) 



Here the density function is separated into a host and an impurity part: cr(A) = ah{A) + 
(l/L)o"jmp(A). Using the dressed energy defined by 

.(A) = f |()(2A - 2) - .1 - /; dA' |^-i£^.(A'), (102) 



we can express the total energy as [^,|4 
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oo 



The dressed charge is determined by the following equation 

which gives ^s(-B) = reflecting SU(2) symmetry. The variation AB is expressed in 

terms of ASh, i.e. AS'/^ = —ah{B)^s{B)AB. As in the preceding sections, we then obtain 
the flnite-size spectrum of the spin sector 

^E-l + j-^E-^ = ^-^{ASnf + ^(LA..)^ (105) 

where we have suppressed the particle-hole excitation just for convenience and the velocity 
of free fermion v = 1. The velocity Vimp of spin excitation for the impurity part is written 
in terms of the dressed energy 

-1 de{A) 



^""^ (^imp{B) dB 



A=B 



(106) 



One can see from Eqs. (p9|) and ( |105| ) that the flnite-size spectrum of the s-d exchange 



model is essentially identical to that of the Anderson model under 6p = tt/2. Our results 
are in agreement with those in [0] for the s-d exchange model. 



V. SUMMARY 

We have studied the flnite-size corrections to the energy spectrum for the Kondo problem 
based on the Anderson model, the degenerate Anderson model and the s-d exchange model. 
All these models are known to exhibit the local Fermi-liquid properties at low energies. The 
basic picture of local Fermi liquid was already established and conflrmed by the exact evalu- 
ation of various static quantities. A recent CFT approach to the Kondo problem by Affleck 
and Ludwig has clarifled how the Fermi-liquid picture emerges and is described in terms of 
boundary CFT. Our calculations demonstrate that the flnite-size spectrum consistent with 
boundary CFT is derived directly from the microscopic models. We have also pointed out 
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that the critical behavior related to the X-ray absorption singularities will also be described 
by boundary CFT. A natural extension of the present work is to analyze the multichannel 
Kondo problem including the overscreened case. It is in principle possible to work out the 
finite-size spectrum using the Bethe-ansatz solution. We hope to turn to this issue in the 
near future. 
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APPENDIX A: FINITE-SIZE CORRECTIONS TO THE GROUND-STATE 
ENERGY OF THE ANDERSON MODEL 



We present the details of computations of the finite-size corrections to the ground-state 
energy of the Anderson model in Sec. I. Applying the Euler-Maclaurin formula to Eqs.(^ 
and one obtains the finite-size form of the Bethe-ansatz integral equations. The result 
reads 



pM = ^ + yAik) + B\k) r dAa.iBik) - A)a,.(A) 

1 ( B' {k)d-^{B{k) - B' {k)a[{B{k) - A- 



24L2 I af,{A+) cr^(A-) 

1 ra2(A-A+) a2(A-A") 



(Al) 



poo ^ 

..(A) + /^_dA'..(A-A'K(A') + — I ^^^^^^ ^^^^ 

J-oo 24L^ I ph{k+) Ph{k ) 

= A{A) + -^Z(A), (A2) 

with k^ = f3 and A~ = a. We let k~ — oo as well as A"*" — +00 at the end of the 
calculation. The solutions of Eqs. ( |A1|) and ([A2| ) are written as 



f'Uk) = p,m + + ^{i^ + 41^ + 4^ + (A3) 

L 2AL^[ph{k+) ph{k ) ah{A ) ah{A+) } 

rr(\^ .r(^^M^.r ( ^^ M ^ r ^i+(A) a^{A) a, (A) or+(A) l .... 

aL{A) = (Th{A) + -airnp[A) + -— <^ , X + (r. + 7TTT + /A.N f- (A4) 

L 24L^ lp/i(A;+) p/i(A; ) crft(A ) a/i(A+) J 

Here pf{k) and erf (A) (z = 1,2) satisfy 

POO 

pHk) =f^^ + B (k) / dAa^{B{k) - A)a^ (A), (A5) 
a,f(A)+/ dAa2(A-A)a,±(A)+ / dA:ai(A - 5(/c))p,±(A:) = (A6) 

J a J —00 



where p?^ = 0, p^^ = ±B' {k)a[{B{k) - A±), a?^ = ±fi'(A;±)a;(5(A;±) - A) and a. 
±a;(A±-A). 



o± 
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Let us next apply the Euler-Maclaurin formula to Eq.(^. We then obtain the finite-size 
form of the ground-state energy as follows: 

Eo= l{ dkip{k)kpL{k) + 2 r dAx{AMx{A))aL{A)} (A7) 

1 2a;' (A- ) 

"24Lp;,(A;+) + 24Lor^(A-)- 

To derive this we have used the assumed property (p{k~), (p' {k~) — > as A;~ ^ — oo and 
V9(x(A+)), v9'(x(A+)) as A+ — >• -|-oo. Using Eqs.(|A3|) and ( |A4| ) we express the ground- 
state energy in the form 

Eo = ^^o+ 24La,(A-) {^'^'^'^'^ + J^^dk^{k)kP2 W + 2 dA^ixiA))xiA)a, (A)^ 
+ 24Lal(A+) {/-l^^^^^^^^^^^^ + ^ /JdA^(x(A))a;(A)a+(A)} 
+ 2ALp,{k+) {~^ + f_Jkmkpt{k) + 2 |JdA^(x(A)MA)a+(A)} 
+ ^iZ;^(F){Xld^^(^)^'«r(fc) + 2 /JdAy,(x(A))a:(A)ar(A)}, (A8) 



where 



//3 poo 
dk(p{k)kp{k) + 2 dAx{A)ip{x{A))a{A). (A9) 
-oo J a 

It is easily seen from Eqs.(|X5|) and ( |A6| ) that the second and fourth lines of ( [A8| ) vanish in 
the limit k" — > — oo, A+ +oo. Furthermore the first and third lines can be expressed in 
terms of the dressed energy functions (^) and 

1 dssik) 



En = LEq + ^ 



2ALahia) dA 
Here we note the relations 



+ 



A=a 2ALph{(3) dk 



k=l3 



(AlO) 



1 de,{A) 



ah{a) OA 



= -27r, (All) 

k=l3 



A=Q ph{,l3) dk 

which can be checked from Eqs.(P), ([T0|), ([T8| ) and ([T9|). Thus we arrive at Eq.(P) in the 
text. 
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